Composition, structure, and stability of the rutile TiO2(110) surface: oxygen 
depletion, hydroxylation, hydrogen migration and water adsorption 
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A comprehensive phase diagram of lowest-energy structures and compositions of the rutile 
TiO2(110) surface in equilibrium with a surrounding gas phase at finite temperatures and pres- 
sures has been determined using density functional theory in combination with a thermodynamic 
formalism. The exchange of oxygen, hydrogen, and water molecules with the gas phase is considered. 
Particular attention is given to the convergence of all calculations with respect to lateral system 
size and slab thickness. In addition, the reliability of semilocal density functionals to describing the 
energetics of the reduced surfaces is critically evaluated. For ambient conditions the surface is found 
to be fully covered by molecularly adsorbed water. At low coverages, in the limit of single, isolated 
water molecules, molecular and dissociative adsorption become energetically degenerate. Oxygen 
vacancies form in strongly reducing, oxygen-poor environments. However, already at slightly more 
moderate conditions it is shown that removing full Ti02 units from the surface is thermodynam- 
ically preferred. In agreement with recent experimental observations it is furthermore confirmed 
that even under extremely hydrogen-rich environments the surface cannot be fully hydroxylated, 
but only a maximum coverage with hydrogen of about 0.6-0.7 monolayer can be reached. Finally, 
calculations of migration paths strongly suggest that hydrogen prefers to diffuse into the bulk over 
desorbing from the surface into the gas phase. 

PACS numbers: 68.43.Fg, 68. 43. Be, 68.47.Gh, 82.65. -l-r 



I. INTRODUCTION 

Over the past years, the (110) surface of Ti02 in the ru- 
tile structure has become one of the most popular model 
systems for fundamental surface science studies of transi- 
tion metal oxides^i It is the thermodynamicaUy most sta- 
ble crystal face of Ti02 and therefore represents the dom- 
inating facet of rutile crystallites.^ Stoichiometric, single 
crystal (1x1) surfaces may be easily prepared, and most 
experimental surface science techniques can be applied 
without difhcultiesji The interest in TiO2(110) surfaces 
is furthermore driven by many technological applications 
of Ti02, ranging from pigments, coatings, electronic de- 
vices, implants, gas sensors, photochemical reactions to 
catalysis. In all of them, the surface properties of Ti02 
play a crucial rolei 

One of the most important properties of Ti02 is that 
it can be easily reduced. The reducibility is essential 
for many applications of Ti02 in heterogeneous cataly- 
sis. Oxide-supported metal-based catalysts with Ti02 as 
part of the support often show a so-called strong metal- 
support interaction (SMSI).'^ Here, the catalytic proper- 
ties of the supported metal clusters are profoundly mod- 
ified by an incorporation of partially reduced TiO^ into 
the boundary areas of the metal particles.^ In many cases, 
the TiOa^-decorated metal particles exhibit a much higher 
catalytic activity for hydrogenation reactions than the 
pure metal itself^*' 

One way to reduce the TiO2(110) surface is to re- 
move surface O atoms. By this process, formally two 
neighboring Ti^+ ions of a vacancy change to a Ti'^^ 



oxidation state. In ultra-high vacuum (UHV) experi- 
ments O vacancies are easily created either by electron 
bombardment, sputtering, or simply by annealing. The 
presence of O vacancies strongly increases the reactiv- 
ity of the surface. Among the many investigated pro- 
cesses, probably the best studied surface reaction is the 
dissociation of water, which has been shown to occur 
at O vacancies (^i^i^ii^ whereas on the stoichiometric, 
well-annealed parts of the surface the water molecules 
stay mostly intact i^^d^d^ ^pi^ere is, however, a limit to 
what extend the surface can be reduced. Oxygen va- 
cancy concentrations are typically in the order of several 
percent,i-'^^ but it is not possible to remove all surface O 
atoms, in contrast to, for example, the rutile SnO2(101) 
surface li^iii 

Alternatively, the TiO2(110) surface can be reduced by 
hydroxylation of the surface O atoms via adsorption of 
hydrogen. While on unreducible oxides, such as MgO, 
a heterolytic dissociative adsorption of H2 with H+ ad- 
sorbing on 0^~ and H~ on the metal cations is preferred, 
it is more favorable for reducible oxides, such as Ti02, to 
form only 0H~ groups. In the latter case the excess elec- 
trons are transferred to the cations, thus reducing Ti'*+ 
to Ti3+. 

The interaction of hydrogen with TiO2(110), how- 
ever, has been much less intensively investigated by 
surface science studies than the reduction by O deple- 
tion. This is rather surprising in view of the impor- 
tance of Ti02 as catalyst component for hydrogena- 
tion reactions and the prospective application of Ti02 
as photocatalyst for the decomposition of water. It 
has been shown that molecular hydrogen does not in- 
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teract strongly with TiO2(110),"'^*'"'^^ while atomic hydro- 
gen readily sticks to the surface O atoms. No 
Ti-H vibrations could be detected with high resolution 
electron energy loss spectroscopy (HREELS).^^ Interest- 
ingly, also the reducibility of TiO2(110) with hydrogen 
is limited. Recent scanning tunneling microscopy (STM) 
measurements^^, revealed that even after very high ex- 
posures to atomic hydrogen the surface cannot be satu- 
rated. Only a maximum coverage of about 0.7 monolay- 
ers could be achieved. The scattering of thermal energy 
He atoms (HAS) showed that the hydroxyl groups do not 
form an ordered overlayer.^^ Performing a thermal des- 
orption spectroscopy (TDS) experiment by monitoring 
the He atom reflectivity of the surface while increasing 
the temperature (He-TDS), two distinct changes in the 
reflectivity at 388 K and 626 K were observed, which in- 
dicate structural rearrangements of the surface.— These 
rearrangements must involve the loss of the hydrogen 
atoms, since afterwards the surface was H-free, as seen 
by HAS and HREELS. However, in conventional TDS 
no desorption of H2, and only a very small amount of 
H2O was detected)^ even after heating the surface up to 
650 K. This led to the conclusion that the H atoms rather 
diffuse into the bulk than desorb from the surface 

As for experiment, also theoretical studies of 
reduced TiO2(110) have mainly focused on O de- 
flcient surfaces. While the properties of O va- 
cancies on TiO2(110) have been investigated 
extensively23^i25^i27^i29^^^^^^ (g^g Ref.Sfor 

a recent review) , we are aware of only two recent studies 
addressing specifically the adsorption of hydrogen.i^^i^ 
The determination of binding energies, adsorption sites 
and some reaction barriers, however, was restricted 
in both cases to slabs with a (1x1) periodicity, thus 
assuming H coverages of one and more H atoms per 
surface unit cell. 

Finally, next to the formation of O vacancies and the 
hydroxylation of the surface by hydrogen adsorption, also 
the interaction with water has to be taken into account. 
The adsorption of water, either molecular or dissocia- 
tive, does not lead to a reduction of the surface. But 
since water is always present, even in well-controlled 
UHV experiments, one has to account for residual hy- 
droxyl groups which affect other adsorption and reaction 
processes.- Water adsorption on TiO2(110) has been in- 
vestigated extensively both experimentally and theoreti- 
cally. Among the experimental studies the common view 
is that at all coverages water adsorbs molecularly on the 
ideal terraces of the TiO2(110) surface and only disso- 
ciates at defectsiiiiiiiS'^^ However, it should be noted 
that in all these experiments always a small amount 
of dissociated water was present as seen by XPS and 
HREELS and indicated by a high-temperature tail in 
TDS. These signatures were naturally attributed to wa- 
ter molecules dissociated at O vacancies, but it would 
be difficult to distinguish them from a situation in which 
water initially would adsorb dissociatively at very low 
coverages and molecularly afterwards as suggested by 




FIG. 1: Atomic structure of the stoichiometric rutile 
TiO2(110) surface. Oxygen atoms are shown in red, Ti atoms 
in gray. The two- and threefold-coordinated O sites and the 
fivefold-coordinated Ti sites are labeled by 0(2), 0(3) and 
Ti(5), respectively. 



some early studies In contrast to the consensus 
among the experimental studies the results from theo- 
retical calculations are very contradictorily. Early inves- 
tigations predicted water dissociation at all coverages in 
complete disagreement with experiment. Only in some 
more recent calculations partial dissociated or molecu- 
lar structures of water were found to be lower in en- 
ergy. Overall, the theoretical studies are almost evenly 
divided whether molecular, dissociative or a partial dis- 
sociative adsorption of water (strongly depending on the 
coverage) is predicted as the energetically most stable 
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The aim of the present paper is to explore the com- 
petition between the reduction of the TiO2(110) surface 
via O vacancy formation or hydrogen adsorption and the 
non-reductive interaction of the surface with water in a 
comprehensive and systematic way. For a series of sur- 
face models with different O defect configurations and for 
a wide range of both hydrogen and water coverages we 
determine the total energies and the fully relaxed atomic 
structures using first-principles density functional theory 
(DFT). In order to extend the zero temperature and zero 
pressure DFT results to relevant environmental situa- 
tions, such as UHV or reaction conditions in heteroge- 
nous catalysis, and to identify the thermodynamically 
most stable surface structures and compositions depend- 
ing on the experimental conditions, we assume that the 
surfaces are in thermodynamic equilibrium with a sur- 
rounding gas phase at a given temperature T and finite 
partial pressures p. To account for the exchange of oxy- 
gen and hydrogen between the surface and the gas phase, 
appropriate chemical potentials fiQ(T,p) and iiYi{T,p) 
are introduced4^i^i^i^ By minimizing the Gibbs free 
surface energy as function of the chemical potentials, sur- 



face phase diagrams of the most stable surface structure 
and composition are constructed depending on both tem- 
perature and partial pressures j^i^i^i^ 

Within this context the main focus of the present study 
will be on the hydroxylated surface. In addition to the 
thermodynamic considerations we investigate the kinetic 
behavior of hydrogen atoms and water molecules on the 
TiO2(110) surface. Energy barriers for various surface 
processes related to the migration and desorption of hy- 
drogen and to the dissociation of water are calculated. 
Our results support the recently proposed suggestion 
that hydrogen atoms, instead of being desorbed from the 
surface at higher temperatures, migrate into the bulk. In 
addition, we find a rather small barrier for water dissocia- 
tion. However, in contrast to some previous studies, our 
results indicate that the molecular adsorption of water 
is preferred over dissociation in the monolayer coverage 
limit. 



II. COMPUTATIONAL APPROACH 
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FIG. 2: Convergence of the unrelaxed (upper panel) and re- 
laxed (lower panel) TiO2(110) surface energy with slab thick- 
ness; note the different energy scale. The solid black line in 
the lower panel represents the result from fully relaxed slabs. 
The calculation with the atoms in the bottom two trilayers 
fixed at the bulk positions is shown by the red dashed line. 



The DFT calculations for the different TiO2(110) 
surface structures, as well as the bulk and molecu- 
lar reference energies, have been carried out using 
the Car-Parrinello Molecular Dynamics (CPMD) 
code.^^'^^ The gradient-corrected Perdew-Burke- 
Ernzerhof functional^ PBE was used to describe the 
exchange and correlation effects. Vanderbilt ultrasoft 
pseudopotentials^ were employed together with a plane 
wave basis set to represent the Kohn-Sham orbitals. A 
plane wave cut-off energy of 25Ryd was sufhcient to 
get well converged results for structures and energetics. 
In particular, increasing the cut-off to 30Ryd changed 
the adsorption energies for water molecules by less 
than 0.02 eV. The Ti pseudopotential was constructed 
from an ionic 4s^ configuration, and the 3s and 3p 
semicore electrons were treated as full valence states. 
Since very large supercells were used (see below), the 
/c~point sampling was restricted to the F-point. Spin 
polarization was included for all systems with an odd 
number of electrons and for all calculations with oxygen 
vacancies. All configurations were relaxed by minimizing 
the atomic forces. Convergence was assumed when the 
maximum component of the residual forces on the ions 
was less than 0.01 eV/A. 

With this computational setup we find for the opti- 
mized Ti02 bulk lattice parameters values of a = 4.649 A, 
c = 2.966 A and u = 0.305, which compares very well to 
previous GGA calculation924i2ii32^^^^57 ^^^^ g^- 

periment (a = 4.594 A, c = 2.959 A and u = 0.305, see 

For the calculation of the energetics of O vacancy for- 
mation as well as hydrogen and water adsorption, the to- 
tal energies of the isolated O2, H2 and H2O molecules are 
needed. While H2 and H2O are reasonably well described 
within PBE/DFT, it is a well-known deficiency of all lo- 
cal and semilocal functionals (such as PBE) that they 
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FIG. 3: Convergence of the hydrogen adsorption energy i?ad 
with slab thickness at full monolayer coverage. The atoms in 
the bottom two trilayers were fixed at the bulk positions. The 
solid black and the dashed red lines distinguish calculations 
without and with saturation of the broken surface bonds at 
the bottom of the slabs with pseudo atoms of nucleus charge 
of +4/3 and +2/3, respectively. 



strongly overbind the O2 molecule i^^i^^'^°i'^^'^°i^'^ In our 
setup we find an O2 binding energy of 5.87 eV, which is 
about 0.4 eV lower than the fully converged PBE value of 
6.24 eV,^"*'^^ which is in turn about 1 eV larger than the 
experimental valued of 5.26 eV (with zero point vibra- 
tional energy removed in the harmonic approximation). 

In order to circumvent errors introduced by a poor 
description of the O2 molecule, we do not employ the to- 
tal energy of the O2 molecule from the DFT calculation 
as reference, but deduce it via a thermodynamic cycle 
from the total energies of the H2 and H2O molecules to- 
gether with the experimental value of 2.51 eV— for the 
formation energy of water from H2 and O2 (also taking 
into account the corrections due to zero point vibrations) . 
With this reference energy for the gas phase O2 molecule 
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FIG. 4: Local density of states (LDOS) of Ti and O atoms 
in bulk Ti02 (black solid lines), of the fivefold- and twofold- 
coordinated Ti(5) and 0(2) surface atoms at the bottom of a 
four trilayer TiO2(110) slab (top two layers relaxed, bottom 
two layers fixed at bulk positions, red dashed lines), and after 
saturation of the bottom Ti(5) and 0(2) surface atoms with 
pseudo- hydrogen atoms (see text, blue dashed-dotted lines). 



we obtain a formation energy of bulk Ti02 from metallic 
bulk Ti and O2 of 9.60 eV, which is in excellent agreement 
with the experimental formation enthalpy of 9.73 eV.— 

The second problem we have to address is how reliable 
can we expect the DFT calculations to be when describ- 
ing the reduced state of an insulator/semiconductori^ It 
is a well-known shortcoming of LDA and all commonly 
used GGA functionals that band gaps are substantially 
underestimated due to the insufficient cancellation of the 
self-interaction energy. With our setup we find a band 
gap for bulk Ti02 of 1.86 eV, which is in very good agree- 
ment with previous DFT calculationsri^iiS^i^iS^ but is 
much smaller than the experimental value of 3.03 eV^^ 
as expected. Upon reduction, defect states are created, 
which typically appear in the band gap region. In DFT 
calculations, however, due to the underestimation of the 
band gap, they might be placed erroneously into the con- 
duction band. In such a case, a delocalized conduction 
band state and not the localized defect state would be 
occupied with electrons. In particular, this is partially 
true for O vacancies on TiO2(110). Experiment indicates 
that the O vacancies are responsible for a defect state 
about 0.7-0.9 eV below the conduction band edgcj ^^^'^^i^^ 
A careful survey of the available literature datai^ showed, 
however, that in well converged GGA/DFT calculations 
this defect level is pinned at the bottom of the conduction 
band, giving rise to a defect state which is too delocal- 
ized. 

The problem of both, the band gap and the position 
of the defect levels, could be cured to a large extend by 
employing hybrid density functionalsj ^^i'^^'^^i^^ Unfortu- 
nately, for extensive studies using large supercells as re- 
quired for the present purpose, they are computationally 
still very demanding even if localized basis sets are used. 
On the other hand, the review Ref. [s^ surprisingly re- 
veals that despite the problem of the correct position of 
the defect level, GGA/DFT calculations give quite good 
results for the energetics of the defect formation. While 



the best estimate for the O vacancy formation energy 
using a cluster model and the B3LYP hybrid functional 
is 2.7 eV (which relies on extrapolations to account for 
the finite cluster size and the limited atomic relaxation 
in this calculation, thus adding some uncertainty to this 
value)," the converged DFT result (using the generalized 
gradient functionals PW91, PBE or RPBE) is 3.0 eV-^ 
It is in particular noticeable that despite the underesti- 
mation of the band gap and the too strong delocalization 
of the defect state, the DFT calculation overestimates the 
vacancy formation energy and thus underestimate the re- 
ducibility of Ti02. 

In recent experiments it was shown that the dissocia- 
tive adsorption of water at O vacancies (thereby form- 
ing two neighboring OH groups on the surface) does not 
change significantly the O vacancy induced defect state 
in the band gap (which was seen 0.9 eV below the conduc- 
tion band edge).^*' Also for this hydrogen reduced state of 
the TiO2(110) surface (which can be viewed as adsorbing 
an H2 molecule on the defect-free surface instead of dis- 
sociating a water molecules at an O vacancy) , Di Valenti 
et al.^^ found that the defect state is not separated from 
the conduction band in GGA/DFT calculations. Never- 
theless, despite this problem of accurately describing the 
electronic structure, we will argue in the present paper 
that the GGA/DFT calculations still give quite reliable 
results for the energetics of reduced surfaces via hydrox- 
ylation, similar as we have seen for the surface reduction 
via O vacancy formation. In particular in the regime of 
higher H coverages, when the defect states start to in- 
teract and a defect band is formed which will gradually 
hybridize with the conduction band, the DFT calcula- 
tions will be less and less hampered by the band gap 
problem. 

Figure [1] shows the atomic structure of the stoichio- 
metric TiO2(110) surface. In our calculations all sur- 
face structures were modeled by periodically repeated 
slabs. Slabs for the stoichiometric (110) surface are 
built by a stacking sequence of trilayers with a composi- 
tion of 0-Ti202-0. In the bulk, the Ti and O ions are 
sixfold- and threefold-coordinated, respectively. At the 
surface this coordination is reduced and the characteris- 
tic feature of the (110) surface is the presence of rows of 
twofold-coordinated bridging O ions parallel to fivefold- 
coordinated Ti atom along the [001] direction. Using 
formal ionic charges of +A and —2 for the Ti and O ions, 
respectively, the stoichiometric (110) surface would be 
charge neutral. From a more covalent point of view, the 
(110) surface, as shown in Fig. [1] is created by a cleavage 
of the crystal in which the minimum possible number of 
bonds has been broken. 

The surface energy for the unrelaxed slabs converges 
very fast with slab thickness as shown in Fig. O After a 
full relaxation, however, the surface energy shows strong 
odd-even oscillations with the number of trilayers in the 
slab (see Fig. [5]). A very large number of trilayers would 
be needed to get well converged results. These oscil- 
lations have been noted previously^ i'^°'^^i^^'^^i^^ and are 
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also present in plots of vacancy formation22i2ii^ and wa- 
ter adsorption energies4^i^i^ The odd-even oscillations 
arise from the different symmetry of the slabs which leads 
to a significant change in the hybridization of the O 2p 
and Ti id states in the direction of the surface normaL— 
In slabs with an even number of trilayers the 0~Ti hy- 
bridization causes a modulation of the coupling between 
the trilayers in such a way that a sequence of pairs of 
trilayers is formed in which the coupling of the trilayers 
within the pair is stronger than the interaction between 
pairs;^ For slabs with an odd number of trilayers, in con- 
trast, the atomic relaxations, which are associated with 
this pair formation, are suppressed by the central Ti2 02 
mirror plane and the 0-Ti hybrid orbitals become more 
delocalized over the whole slab.^- The surface energy of 
the slabs with an odd number is higher and the conver- 
gence to the limit of infinite slab thickness is slower than 
when using an even number of trilayers. The odd-even 
oscillations can be largely reduced if not the full slab is 
allowed to relax. In Fig.[5]the surface energy is shown for 
a calculation in which the atoms in the bottom two tri- 
layers were held fixed at their bulk positions. Using this 
setup, quite well converged results are already obtained 
with only four trilayers. 

Based on this observation Thompson and Lewis^ pro- 
posed to use such a slab setup of four trilayers with the 
bottom two layers fixed for any large scale calculation. 
However, the surface properties of a stoichiometric slab 
might be less influenced by the slab thickness than when 
it comes to reduced surface structures. States above the 
valence band (either defect states in the band gap or 
states from the lower edge of the conduction band) are 
more delocalized than the valence states and might be 
more sensitive to the truncation at the bottom of the 
slab. For example, Leconte et alM- have reported that 
for fully hydroxylated slabs with a thickness of two tri- 
layers only 63% of the spin density is localized at the 
hydroxylated top surface layer, but 37% at the bottom 
of the slab. The fully hydroxylated surface is the highest 
reduced state which we will consider in our calculations 
(though, as we will show in Sec. IIIIBl this configuration 
is thermodynamically not stable). Even after fixing the 
atomic positions of the bottom two trilayers to the bulk 
positions, the hydrogen adsorption energy at monolayer 
coverage still shows a slow convergence and noticeable 
oscillations with slab thickness, as can be seen in Fig. [31 

Though Ti02 is a strong ionic insulator, the chemical 
bond has significant covalent contributions. It is there- 
fore worth to think about the Ti02 surface in terms of 
broken surface bonds as one would do for covalent semi- 
conductor surfaces fZi Ti and O have four and six va- 
lence electrons, respectively. To fulfill the octet rule, each 
O atom receives 2/3 electrons from its three nearest Ti 
neighbors and Ti contributes 4/6 electrons to each of its 
six nearest-neighbor bonds. At the TiO2(110) surface the 
fivefold Ti and the twofold bridging O atoms have lost one 
of their nearest-neighbor atoms. Thus, dangling bonds 
are created which are occupied with 2/3 electrons for the 



fivefold Ti(5) and 4/3 electrons for the twofold 0(2). At 
the stoichiometric surface the partial occupation of these 
dangling bonds is removed by autocompensationj^ Al- 
ternatively, we saturate these broken bonds at the bot- 
tom of the slab by introducing artificial atoms with a 
nucleus charge of +4/3 and +2/3 which we place next 
to the Ti(5) and 0(2) atoms, respectively, thus creating 
a more bulk-like environment for these surface atoms^^^ 
The distance between the pseudo-hydrogen and the sur- 
face atoms was determined by a geometry optimization 
in which all Ti and O atoms were kept fixed at the Ti02 
bulk positions. The effect of saturating the broken sur- 
face bonds with the pseudo-hydrogen atoms on the elec- 
tronic structure of the surfaces is illustrated in Fig. H) 
Without saturation, the local density of states (LDOS) 
of the fixed Ti(5) and 0(2) surface atoms (red dashed 
lines) deviates strongly from the bulk behavior (black 
solid lines). The surface band gap is reduced and strong 
changes and shifts in the LDOS peaks due to the auto- 
compensation can be seen. After saturation of the bro- 
ken surface bonds, however, the LDOS of the surface 
atoms (blue dashed-dotted lines) is almost indistinguish- 
able from the bulk. With this saturation of the bottom of 
our slabs we find now also for reduced surface structures 
a much faster convergence of surface properties with slab 
thickness. The oscillations in the hydrogen adsorption 
energy in Fig. [3] have been reduced from 0.1 eV to about 
0.02 eV per adsorbed hydrogen atom. 

Based on these convergence tests we decided to use for 
all further calculations slabs with a thickness of four tri- 
layers, including a saturation of the broken surface bonds 
at the bottom. The upper two trilayers were always fully 
relaxed while the atoms in the lower two layers were fixed 
at the bulk positions. A large (4x2) surface unit cell2^ 
was used in order to be able to study adequately the 
coverage dependence of hydrogen and water adsorption. 
The slab for the stoichiometric surface thus contained 
208 atoms, including the artificial atoms for saturating 
the broken surface bonds at the bottom of the slab. The 
calculated equilibrium bulk values were taken for the lat- 
tice constants parallel to the surface. The slabs were 
separated by a vacuum region of about 13 A thickness, 
which corresponds to the thickness of the slab itself. The 
surface relaxations of the stoichiometric surface are essen- 
tially the same as described by Thompson and Lewis^. 
As shown in Fig. [21 the relaxation energy of the surface 
is quite substantial. The unrelaxed surface energy is re- 
duced by atomic relaxations from 1.64 eV (1.35 J/m^) by 
1.21 eV to 0.43 eV per surface unit cell (0.36J/m_2). For 
comparison, the relaxation energy for the ZnO(lOlO) sur- 
face amounts to only 0.37 eV per surface unit cell.^^ As 
we will see later on, these large surface relaxations have a 
strong impact on the coverage dependence of adsorption 
energies. 

In order to analyze the thermodynamic stability of our 
different surface structures we assume that the surfaces 
can exchange O and H atoms with a surrounding gas 
phase. Assuming thermodynamic equilibrium, the most 
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TABLE I: Vacancy formation energies (in eV) for different 
oxygen vacancy configurations on TiO2(110) in the fimit of 
A/^o = (oxidizing, i.e. oxygen rich conditions). 



Vacancy type 


iSv singlet 


i?v triplet 


0(2)-v 


3.07 


3.02 


0(3)-v 


4.05 


4.00 


Ti02-V 


1.03 





stable surface composition at a given temperature T and 
pressure p is given by the minimum of the Gibbs free sur- 
face energy 7(T, Since we are only interested 
in the relative stabilities of surface structures, we cal- 
culate directly the difference A7(T, p) of the Gibbs free 
surface energy of the defective or adsorbate-covered and 
the stoichiometric, ideal surface according to 

A7(r,p) = i(G:-^(T,p,AiVo,A7VH)-G^L'bmp) 
+ A7VoMo(T,p)-A7VHMH(r,p)) , (1) 

where G|fab ^'^d ^siab ^^'^ the Gibbs free energies of the 
modified and the stoichiometric reference surface config- 
urations, respectively. A is the surface area, ANq, AA^h 
are the differences in the numbers of O and H atoms be- 
tween the two surfaces, and ^o{T,p), ^uiT-^p) are chem- 
ical potentials representing the Gibbs free energy of the 
gas phase with which the O and H atoms are exchanged. 
According to this definition, A7 is negative if the mod- 
ified surface is thermodynamically more stable than the 
stoichiometric surface and positive otherwise. Assum- 
ing that all differences in entropy and volume contribu- 
tions in A7 are negligibler^i^ we approximate the Gibbs 
free energies G^^^l and C^l^ by their respective total en- 
ergies of our DFT slab calculations as usua L^^'^^ Up- 
per bounds for the chemical potentials /j,o and /j,h are 
given by the total energies of their most stable elemental 
phases?^ that is, molecular oxygen \E^^^ and molecu- 
lar hydrogen \E^^^\, respectively. These upper bounds 
are taken as new zero point of energy by introducing 

A^o = ^lo- l^^moi and A^h = m~ \E^Iv A lower 
bound for A/^o is given by minus half of the formation 
energy of bulk TiO^, i.e. eJ'''^ = E^,^ + E^l, - i?™^ 
(here E^^ and E'^^^^ are the energies of one bulk unit 
cell of Ti02 and metallic Ti, respectively) ^Si^ for which 
we have taken the theoretical value of 4.80 eV from our 
PBE/DFT calculations. The chemical potential can be 
related to experimental temperature and pressure con- 
ditions by using experimental thermochemical reference 
data or by applying the ideal gas equationi^Si^iSi 

The transition state search for the dissociation, des- 
orption and migration processes was conducted with 
the nudged elastic band (NEB)^^ and the dimer 
method.— i^i^i^^ Fully relaxed configurations were cho- 
sen as initial and final state of the NEB calculations. 
Throughout twelve images were used which were con- 
nected by springs with a fixed spring constant of about 



20eV/A2. The two images in the dimer calculations were 
separated by 0.01 A in configuration space and the trial 
steps for translation and rotation were 0.01 A and 10°, re- 
spectively. In all calculations first a good approximation 
of the transition path was determined with NEB which 
was then refined with a dimer method run for a precise 
location of the transition state. 



III. RESULTS AND DISCUSSION 
A. Oxygen Vacancies on TiO2(110) 

We have investigated three different types of oxygen- 
related defects. Vacancies are created by removing a sin- 
gle bridging oxygen 0(2) atom, a threefold-coordinated 
0(3) atom from the surface layer, or a full Ti02 unit of 
two neighboring bridging 0(2) atoms and the underlying 
Ti cation, see Fig.[TJ With our choice of (4x2) supercells 
this corresponds to a vacancy concentration of 1/8 mono- 
layer and a defect separation of about 12 A. The vacancy 
formation energies for the missing 0(2) and 0(3) atoms 
and the Ti02 unit are given by^ 

= <V + \eZI^ - KfL + A/^o (2) 

and 

Ev - -t^slab + ^bulk - ■'^slab : l-^) 

respectively, where E2~f^ , -^siab^"^ ^'^'^ ^sfab ^''^ ^^ic to- 
tal energies for the defective and the stoichiometric slab, 
respectively. The results in Table H] show that the 0(2) 
and 0(3) vacancies are slightly more stable in the triplet 
than in the singlet state, in agreement with previous 
calculations f2^j^i^i2i As one would expect from the lo- 
cal coordination, the formation of an 0(2) defect is much 
more favorable than the 0(3) vacancy. Our value for 
the 0(2) vacancy formation energy of 3.02 eV is basically 
identical to the results of Rasmussen et ali^ (3.03 eV) 
and Oviedo et alM- (3.07 eV, singlet state) who have both 
conducted extensive studies on the convergence of EV'^> 
with defect concentration, slab thickness and the num- 
ber of relaxed surface layers. Also Ganduglia-Pirovano et 
al— came to the conclusion in their survey of the avail- 
able literature data that the converged value for the for- 
mation energy of isolated 0(2) vacancies in GGA/DFT 
calculations (using the PW91, PBE or RPBE functional) 
is about 3.0 eV. This supports that our choice of slab 
thickness, relaxation, bond saturation at the bottom of 
the slab, and size of surface unit cell is indeed appropriate 
to give well converged results. 

In Fig. [5] the vacancy formation energies are plotted 
as function of the oxygen chemical potential A/j,o of the 
environment. Assuming thermodynamic equilibrium of 
the surfaces with an O2 gas phase we have converted the 
chemical potential into a pressure scale for a temperature 
of 800 K and 1200 K using the ideal gas equationi^i^ 
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FIG. 5: Defect formation energy for O vacancies at the 
TiO2(110) surface as function of the oxygen chemical poten- 
tial A/io- In the top x-axis, the chemical potential has been 
translated into a partial pressure scale po^ of molecular oxy- 
gen at 800 K and 1200 K. 



Since we have neglected the entropy contributions to 
the Gibbs free energy of the defective and stoichiometric 
slabs, we cannot capture with our description the for- 
mation of thermal oxygen vacancies which are due to 
the gain in configurational entropy by the random dis- 
tribution of the defects on the surface. Therefore, at 
lower temperatures we will always find the stoichiomet- 
ric, defect-free surface to be the most stable one. We 
can only predict the appearance of structural vacancies 
which are formed if the gain in entropy by bringing oxy- 
gen into the gas phase outweighs the binding energy to 
the surface. Having this in mind we have to expect that 
we will overestimate the temperature at which significant 
amounts of O vacancies (in the percent range) will form. 
From Fig. [5] we see that at UHV conditions (base pressure 
of about 10~^° mbar) the onset for the formation of struc- 
tural 0(2) vacancies is around 1200 K (i.e., the chemical 
potential at which the vacancy formation energy be- 
comes zero) . This is in quite good agreement with exper- 
imental observations fi. taking into account that an uncer- 
tainty of the vacancy formation energy of a few tenths of 
an eV from the DFT calculations translates into an er- 
ror bar for this temperature in the range of 100-200 K. 
Rasmussen et al^ and Oviedo et alM^ have shown that 
the interaction between the 0(2) vacancies is strongly 
repulsive and that the O vacancy formation energy in- 
creases rapidly with higher defect concentrations. Thus, 
structures with higher O defect concentrations will ap- 
pear farther to the left in the phase diagram Fig. [S] so 
that increasingly reducing conditions are needed to reach 
defect concentrations beyond a few percent. 

On the other hand, quite surprisingly, we find that the 
formation energy for the Ti02 defects is much lower than 
for the 0(2) and 0(3) vacancies over a wide range of the 



TABLE II: Adsorption energies of H atoms (in eV/atom) 
on TiO2(110). Configurations with negative values of E]^ are 
energetically unstable towards the desorption of H2 molecules. 
A^H is the number of hydrogen atoms in the (4x2) surface 
unit cell. The arrangements of the H atoms on the surface 
are illustrated in Fig. |6l 
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-0.012 
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O chemical potential. Already at temperatures around 
800 K (at UHV conditions) it becomes thermodynami- 
cally more favorable to convert the 0(2) vacancies into 
Ti02 defects. A similar result was recently obtained for 
the ZnO(lOiO) surface,— where it was found that ZnO 
dimer vacancies are thermodynamically more stable than 
O vacancies for most experimental conditions. While for 
ZnO a strong suppression of O vacancies on the (1010) 
surface can be expected, since Zn has a very low vapor 
pressure and can easily desorb from the surface (indeed, 
desorption of Zn is regularly observed in TDS experi- 
ments of ZnO), it is more difficult in the case of Ti02 
to remove the Ti ions from the surface since the kinet- 
ics of the process will play a much larger role. How- 
ever, one has to consider that the O vacancies, created 
on the TiO2(110) surface in UHV experiments by an- 
nealing, sputtering or electron bombardment, are only 
present due to kinetic limitations, and they will be elim- 
inated under ambient conditions not only by a dissocia- 
tive adsorption of O2 and water, but there is also a strong 
driving force to transform them into Ti02 vacancies, for 
example, by a mechanism in which Ti interstitials are 
createdi^ 

B. Hydrogen Adsorption on TiO2(110) 

In order to investigate the interaction of hydrogen with 
the TiO2(110) surface we have calculated the total ener- 
gies of various hydroxylated TiO2(110) surface structures 
taking into account both, different hydrogen coverages 
and adsorption sites for the H atoms. Our results for the 
hydrogen adsorption energies i?^ per H atom 

El = ^ {e'A + ^E^Z, - i?lT'(^H)) , (4) 

where ^'^^^'''^(^h) is the total energy of the slab calcu- 
lation with iVu adsorbed H atoms, are summarized in 



8 



(a) 0(2)+0(2) 0(2)+0(3) 0(2)+Ti(5) 




1_1 

tl ft I 1 t t 



4H-e 5H 6H-a 6H-b 

t | t | t | l | | t | t | t | I | t | t | t | I | t | t | t | l 
t I t t 1 t I 



FIG. 6: Schematic illustration of some configurations of the 
TiO2(110) surface with adsorbed H atoms, (a) Configurations 
with two adsorbed H atoms. Ti, O and H atoms are depicted 
as small gray, large red and small white spheres, respectively, 
(b) Configurations with higher H coverage. The arrows in- 
dicate the tilt of the OH groups after If adsorption. The 
corresponding adsorption energies iJad are given in Table HIl 



Table [TTl Single H atoms adsorb preferentially ontop of 
the bridging 0(2) atoms, whereas the 0(3) and Ti(5) 
sites are significantly higher in energy. The OH groups, 
which form upon H adsorption, do not remain upright 
but break the mirror symmetry of the surface by tilting 
about 20° (isolated OH groups) and up to 50° (hydroxy- 
lation of all bridging O atoms) with respect to the surface 
normal. 

For pairs of H atoms we find that the homolytic ad- 
sorption on two neighboring 0(2) sites (thereby reducing 
the surface) is more stable than the heterolytic adsorp- 
tion on 0(2) and Ti(5) (for the atomic configurations, see 
Fig. [5^). The heterolytic adsorption is even energetically 
unstable towards desorption of H2. This preference of 
the homolytic adsorption of hydrogen is expected for an 
easily reducible oxide. Table HTl also shows that the inter- 
action between neighboring OH groups is repulsive. The 
adsorption energy per H atom decreases from 0.56 eV for 
isolated OH groups to 0.40 eV for neighboring OH pairs. 

By comparing the energy of the surface structure with 
two neighboring OH groups to the structure with an 0(2) 
vacancy, we can directly deduce the energy gain for dis- 
sociatively adsorbing a water molecule at isolated 0(2) 
defects by applying a thermodynamic cycle. For the wa- 
ter adsorption energy at 0(2) vacancies we find a value of 
1.31 eV, which is slightly larger than previously reported 
results (0.94, 0.97 and 1.10 eV according to Refs. T.'s'.'ss', 
respectively), but agrees well with an estimated value of 




0.2 0.4 0.6 0.8 



hydrogen coverage [ML] 

FIG. 7: Hydrogen adsorption energy iJad as function of the 
hydrogen coverage of the surface. 



about 1.4 eV deduced by a Redhead analysis^ of a wa- 
ter desorption peak at 520 K in TDS experiment which 
was assigned to the recombination of neighboring OH 
groups f2£ According to our NEB calculations this pro- 
cess is without barrier— so that the activation energy 
for desorption is given by the adsorption energy. Since 
the previous calculations of Refs. [tIIsIIssI were done with 
(2x2) surface unit cells for which there is still significant 
interaction between O vacancies and pairs of hydroxyl 
groups (which can be seen, for example, from our results 
in Table ITT)) . we attribute our slightly higher value for 
the binding energy of water at O vacancies to our larger 
(4x2) surface unit cell which represent better converged 
results with respect of describing isolated defects and OH 
pairs. It is interesting to see that we underestimate the 
binding energy of water in the O vacancies. Since the 
PBE/DFT calculations tend to overestimate the O va- 
cancy formation energy^^^, this means that also the en- 
ergy of the reduced surface via hydroxylation is slightly 
too high. Thus, as in the case of surface reduction via 
O depletion, also the reducibility of Ti02 via hydrogen 
adsorption is underestimated in the PBE/DFT calcula- 
tions. 

With increasing hydrogen coverage, 0(2) remains the 
preferred adsorption site up to the coverage of a full 
monolayer. No significant amounts of H atoms will be ad- 
sorbed on 0(3) and Ti(5) sites. Neighboring OH groups 
in the [001] direction tend to tilt in opposite directions 
(see Fig. whereas the coupling between [001] rows 
(concerning the tilt direction of the OH groups) is rel- 
atively weak (with a minor preference of the same tilt 
direction of OH groups in [110] rows). The ground state 
configuration of the fully hydroxylated surface, therefore, 
will have a (2x1) symmetry. However, since the barrier 
for flipping the orientation of an OH group is very small, 
they will be fully orientally disordered at room tempera- 
ture. 

Though 0(2) remains to be the preferred adsorption 
site also with increasing H coverage, the adsorption en- 
ergy iJ^ per H atom, i.e. the feasibility to further reduce 
the TiO2(110) surface, decreases strongly upon adsorp- 
tion of more hydrogen (see Fig. [7]). At full monolayer 
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FIG. 8: Surface Gibbs free energy A7 of TiO2(110) surfaces 
with different hydrogen coverages as function of the hydrogen 
chemical potential A/xh- In the top x axis, the chemical po- 
tential has been translated into a partial pressure scale 
of molecular hydrogen at 400 K. Thermodynamically unstable 
surface structures are represented by dashed lines. 



coverage E^^ is still positive. However, if the interac- 
tion energy of hydrogen with the TiO2(110) surface is 
not expressed as binding energy per H atom but as en- 
ergy gain per surface area it becomes obvious from the 
phase diagram Fig. [5] that the fully hydroxylated surface 
is thermodynamically unstable. The Gibbs free surface 
energy A7 for a hydrogen monolayer (this also applies 
for the 0.75 ML coverage) is higher than the surface en- 
ergies of lower coverages over the whole range of the hy- 
drogen chemical potential. I.e., the energy of the fully 
hydroxylated surface can be always lowered by desorbing 
H2 molecules and reducing the H coverage. The highest 
coverage that can be reached in thermodynamic equilib- 
rium in the limit of zero hydrogen chemical potential (i.e. 
at low temperature and high hydrogen partial pressure) 
is around 60%~70%. This is in excellent agreement with 
recent experimental observations where by exposing the 
TiO2(110) surface to atomic hydrogen no H coverage be- 
yond 70% could be obtained.^ The experimental setup 
of exposing the surface to atomic hydrogen is an almost 
ideal realization of the theoretical assumption underlying 
the calculation of the phase diagram Fig. [51 namely that 
the surface is in thermodynamic equilibrium with a sur- 
rounding H2 gas phase: no kinetic limitations are present 
and H atoms can adsorb and desorb without encounter- 
ing any barrier. 

In conclusion, this remarkable agreement between the- 
ory and experiment suggests that despite the problem 
of the band gap and the correct position of the defect 
levels (which we analyzed and discussed in Sec. [TT|, the 
PBE/DFT calculations seem to be quite able to capture 
accurately the subtle balance between the energy cost 
to dissociate H2 molecules versus that of reducing the 
TiO2(110) surface. 



TABLE III: Activation energy _Eact (in eV) for different mi- 
gration paths of H atoms on the hydroxylated TiO2(110) sur- 
face. All calculations are based on a fully hydroxylated sur- 
face with eight H atoms in the (4x2) surface unit cell. The 
temperature Tkct (in K), at which the onset of the processes 
can be expected, has been estimated from the activation en- 



ergy _Eact by applying a standard Redhead analysis 



with 



a frequency factor of 1-10 s and a heating rate of 1 K/s. 



Migration Path 



0(2) ^ 

0(3) - 
Ti(5) - 
0(2) - 
0(2) - 
0(2) - 
2*0(2) 



Osub 

(-.a 

'-'sub 

* Osub 

0(3) 
Ti(5) 
0(2) 

^ H2(gas) 



2.56 
0.93 
1.52 
0.63 
1.54 
1.21 
1.85 



963 ± 197 
361 ± 72 
581 ±117 
248 ± 49 
588 ± 119 
466 ± 93 
703 ± 143 



C. Migration of Hydrogen Atoms 

Upon heating the hydroxylated TiO2(110) surface a 
quite unexpected behavior was observed ;^^i^^ Monitoring 
the He atom reflectivity of the surface revealed two dis- 
tinct maxima at 388 K and 625 K (with a small shoulder 
at 560 K) at which the surface undergoes structural rear- 
rangements. These rearrangements must involve the loss 
of H atoms since the surface is H-free above 650 K. When 
the Redhead formula&S with a frequency prefactor of 
10^^^ s^^ is applied, the two temperatures can be trans- 
lated (with an uncertainty of about 20%) into activation 
barriers of 1.03 eV and 1.68 eV, respectively, for the pro- 
cesses which are involved in the structural changes. In a 
first assignment the two processes were assumed to be the 
desorption of hydrogen from the Ti(5) and 0(2) sitesi^, 
respectively. Based on our results on the hydrogen ad- 
sorption energy for the different surface sites, however, 
this interpretation can be clearly ruled out. Surprisingly, 
in a subsequent conventional TDS measurement, no des- 
orption of H2 at all and only a very small amount of water 
could be detected. This led the authors to the conclu- 
sion that instead of desorption taking place, the H atoms 
migrate into the bulk.^^ 

Aiming at identifying the processes which lead to the 
structural changes of the surface and the loss of the H 
atoms we have calculated the activation barriers for dif- 
ferent migration paths of H atoms on the hydroxylated 
TiO2(110) surface. Some of these migration paths have 
been considered already in a previous DFT study by Yin 
et al.^ in which, however, only the smallest possible 
(1x1) surface unit cell with two adsorbed H atoms was 
used. In our calculations we start from the fully hydrox- 
ylated surface with eight H atoms adsorbed ontop of the 
bridging 0(2) atoms in a (4x2) surface unit cell. In the 
first step we examine how much the energy increases if 
one of the H atoms from the 0(2) sites is placed ontop 
0(3), Ti(5), or subsurface. As we can see from Table HIl 
all three sites are energetically almost degenerate and 
only about 0.3 eV higher in energy than the 0(2) site 
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TABLE IV: Literature survey of calculated water adsorp- 
tion energies E^^ (in eV/molecule) at monolayer coverage for 
molecular ('mol'), mixed dissociative/molecular ('mix') and 
dissociative ('diss') adsorption on TiO2(110). 



Functional 


Reference 


mol 


EZ mix 


EZ diss 


gpgg»4,»!, 


Ref . 38^ 


0.82 




1.08 


PW91-- 


Ref. 40 


1.13 




0.98 




Ref. 41 


0.99 


1.10 


0.91 


RPBE**^ 


Ref. 44 




0.64 


0.37 


RPBE**^ 


Ref. 45 


0.52 


0.53 


0.34 


PW91-'^ 


Ref. 46 


1.09 


1.05 


0.91 


PW91— 


Ref. 47 


1.01 


0.95 


0.90 


pbe'^* 


this work 


0.82 


0.77 


0.63 



(note that the values in Table |TT] are adsorption energies 
per H atom. Therefore, to get the energy difference be- 
tween two structures with different arrangement of the H 
atoms, the E^^ values have to be multiplied by the num- 
ber of adsorbed H atoms, i.e. eight, before subtracting 
the energies) . This energy difference between the adsorp- 
tion sites increases when the H coverage decreases. At 
1/4 monolayer H coverage the 0(3) and Ti(5) sites are 
higher in energy by 0.52 eV and 1.09 eV than the 0(2) 
site, respectively (see Table HI)) . 

Returning back to the fully hydroxylated surface we 
find that the transition barrier for migration of an H 
atom directly from an 0(2) site to a subsurface O atom 
is quite high (see Table IIIip . whereas diffusion to the 
0(3) and Ti(5) sites is much more favorable. On the 
other hand, the activation energies for transferring an H 
atom from 0(3) and Ti(5) to a subsurface O atom are 
moderate. Overall, the energetically most favorable mi- 
gration path is from 0(2) via 0(3) to a subsurface site 
with barriers of only 0.63 eV and 0.93 eV, respectively. 
In contrast, we find for the desorption of H atoms from 
two neighboring 0(2) sites via recombination and forma- 
tion of H2 molecules an activation barrier of 1.85 eV. This 
implies that it is indeed much easier for the H atoms to 
diffuse into the bulk instead of desorbing from the sur- 
face. The activation barrier of 0.93 eV also corresponds 
nicely to the temperature of 388 K at which the onset of 
structural changes of the hydroxylated TiO2(110) surface 
was observed in the HAS experiments.-- Altogether our 
results fully support the conclusions of Yin et al?"^ that 
H atoms from a hydroxylated TiO2(110) surface do not 
desorb but prefer to migrate into the bulk instead. 



D. Water Adsorption on TiO2(110) 

From the experimental point of view there is a wide 
consensus that water adsorbs molecularly on defect-free 
terraces of the TiO2(110) surface .iii^i^ii^ In XPS^i and 
HREELS— measurements only small amounts of disso- 
ciated water could be detected under UHV conditions, 
which may be naturally linked to the presence of resid- 



TABLE V: Literature survey of calculated water adsorption 
energies E^^ (in eV/molecule) of isolated molecules for molec- 
ular ('mol') and dissociative ('diss') adsorption on TiO2(110). 
The size of the surface unit cell used in the respective calcu- 
lation is given in the second column. The (2x2) and c(4x2) 
cells correspond to an effective water coverage of 1/4 mono- 
layer, whereas a (4x2) cell represents a 1/8 monolayer water 
coverage. 



Functional 


unit cell 


Reference 


Sa"d mol 


Sa"d diss 




RPBBSi 


(2x2) 


Ref. 7 


0.56 


-0.23 


-fO.79 


RPBE^ 


(2x2) 


Ref. 45 


0.36 


0.56 


-0.20 


PW91**^ 


(2x2) 


Ref. 46 


0.83 


0.71 


+Q.12 


PBB^ 


c(4x2) 


Ref. 88 


0.76 


0.66 


-hO.lO 


RPBE^ 


(2x2) 


Ref. 


0.66 


0.79 


-0.13 


PBE^"* 


(2x2) 


this study 


0.82 


0.75 


+Q.m 


PBE^'' 


(4x2) 


this study 


0.93 


1.04 


-0.11 



ual defects. In TD S^^'^^i^'^ an intense desorption peak 
is observed at 270 K which could be attributed to the 
desorption of a full monolayer of molecularly adsorbed 
water molecules at Ti^+ sites based on XPS, HREELS 
and work function measurements. Interestingly, the peak 
position shifts to higher temperatures if the initial wa- 
ter coverage O is lowered, indicating a repulsive interac- 
tion between the adsorbate molecules. This is unusual 
since for water molecules one would expect an attrac- 
tive interaction due to the formation of hydrogen bonds, 
but the same behavior is also observed for the anatase 
(101) surface;^ The TDS data could be best described 
by fitting them to a first-order kinetics model. Activa- 
tion energies for desorption of (0.74— 0.09 8) eV-ii and 
(0.73—0.076) eV— were obtained, whereas from a modu- 
lated molecular beam studj*i^ a value of (0.83—0.36 6) eV 
was deduced. 

On the other hand, over the past years theoreti- 
cal studies on the adsorption behavior of water on 
TiO2(110) came to very contradictory conclusions. In 
early GGA/DFT studies a much lower adsorption en- 
ergy for molecular than for dissociative adsorption^i^S 
(at full monolayer coverage) was obtained (it should be 
noted, however, that the binding energy for the molecu- 
lar adsorption was underestimated since only a symmet- 
ric adsorption of the molecules was considered, which is 
not the most stable geometry, see Fig. [H]). Subsequently, 
Bates et al.^'^ found a slightly higher adsorption energy 
for molecular monolayers, whereas Lindan et aZ^ii^i^ 
proposed a mixed molecular/dissociated structure as the 
most stable state. In a Car-Parrinello molecular dynam- 
ics study "'•^ it was observed that water does not dissociate 
on the perfect surface, suggesting that the dissociation is 
hindered by a larger dissociation barrier. Later static 
calculations of the activation energy for water dissoci- 
ation supported this point of view 4^ In a quite recent 
study, however, again molecular adsorption was found to 
be the most stable adsorption mode not only for water 
monolayers but even at lower water coverages 4^ 

The binding energies E'^ from these studies for molec- 



TABLE VI: Water adsorption energies i?"^ (in eV/molecule) 
on TiO2(110) for difTerent coverages taking into account 
molecular(M) and dissociated (D) states of the water 
molecules with downward (d) or upward (u) orientation. TVw 
is the number of water molecules in the (4x2) surface unit cell. 
Parenthesis denote pairs of water molecules as illustrated in 
Fig. [5] The structures with four water molecules consist of 
two water pairs filling one [001] row and leaving every second 
one empty. At full monolayer coverage, both [001] rows are 
equally filled with water pairs. The periodicity of the result- 
ing structure is reported in the 'Config' column. 





Config 




Af„ 


Config 




Ead 


1 


single water 
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double pairs along [001] 




Dd 


1.04 




2(DdD„) 


-(2x2) 


0.67 






u.yo 




2 (DdDd) 


-(1x2) 


0.64 
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0.97 
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0.80 
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0.96 
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0.84 
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-(1x2) 


0.86 
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quadruple pairs 






(DdD„) 


0.94 




4(DdD„) 


-(2x1) 


0.61 




(DdDd) 


0.90 




4 (DdDd) 
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0.63 




(DdM„) 


0.97 




4(DdM„) 


-(2x1) 


0.77 
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0.96 
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-(2x1) 


0.77 
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0.90 
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0.88 




4 (Md Md 


-(1x1) 


0.82 




(Md M„) 


0.88 












(Md Md) 


0.90 











ular, mixed molecular/dissociative, and full dissociative 
adsorption at monolayer coverage are summarized in Ta- 
ble IIVI The situation becomes even more confused if the 
published results for the adsorption of 'isolated' water 
molecules are compared (see Table |Vl . A survey of the 
publications shows that the results for the water bind- 
ing energies depend sensitively on the computational 
setup, in particular the slab thickness and how the re- 
sults are extrapolated to the limit of infinite thick slabs, 
as well as on the GGA functional which was used in the 
calculations (though the result reported in Reff2 on the 
dissociative water adsorption is probably due to an er- 
ror in the calculations since it was never reproduced in 
more recent studies). From embedded cluster calcula- 
tions using wave function based methods (Hartree-Fock, 
B3LYP and MP2) it has been suggested that overall 
GGA/DFT might overestimate the stability of the disso- 
ciated state of the water molecules compared to a molecu- 
lar adsorption42. However, embedded cluster calculations 
have to be taken with considerable caution. As we will 
see below, surface re-relaxations, not only of the nearest 
neighbor atoms of the adsorption site but also including 
atoms in the next surface unit cells, make up a major 
contribution to the binding energy of the molecules, in 
particular for the dissociated state. Due to the limited 
cluster size it is questionable how well such structural 
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FIG. 9: Schematic illustration of water pairs adsorbed on the 
TiO2(110) surface. Ti, O and H atoms are depicted as small 
gray, large red and small white spheres, respectively. The 
labels are defined in Table IVTl 



effects have been captured in this study since the hy- 
droxylated bridging O atoms created by the dissociation 
of a water molecule were already located at the boundary 
of the cluster i-? 

The main intention of our study is not to resolve the 
controversy between the different DFT studies on the 
molecular or dissociative nature of the adsorption of wa- 
ter, but to evaluate the relative stability of surface hy- 
droxyl groups created by the non-reductive dissociation 
of water compared to both molecular water adsorption 
and to hydroxylation of the surface by reductive hydro- 
gen adsorption which we will represent in terms of a 
phase diagram. The general features of the phase dia- 
gram and the qualitative position of the phase bound- 
aries only depend on the overall interaction strength of 
the water molecules with the TiO2(110) surface but not, 
in cases where the energy difference between molecular 
and dissociative adsorption is small as for single, isolated 
water molecules, whether the molecular or the dissociate 
state turns out to be more stable. However, since our 
calculational setup with a large surface unit cell and a 
passivated slab seems to be quite reliable to give well con- 
verged results, we explored the question whether water 
dissociates on the ideal terraces of the TiO2(110) surface 
or not a little bit further than it would have been neces- 
sary for the construction of our surface phase diagram. 

The most favorable adsorption site for water molecules 
and hydroxyl groups is ontop of the fivefold-coordinated 
Ti(5) cations, thereby restoring their sixfold-coordination 
as in the bulk. Water molecules lie relatively flat on 
the surface and form a weak hydrogen bond to one of 
the bridging 0(2) atoms. Hydroxyl groups take a more 
upright position and orient themselves along a Ti(5)- 
0(3) bond, as long as they cannot form a hydrogen bond 
to an adsorbed molecule on a neighboring Ti(5) site (see 
Fig. [9]). 

For single water molecules in a (4x2) surface unit cell 
we find that dissociative adsorption is preferred com- 
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pared to the molecular state by 0.11 eV, as can be seen 
from the compiled results in Table I VII This is in con- 
trast to the results of Ref. but in agreement with 
Ref. '^'^,'45'. On the other hand, at full monolayer cov- 
erage molecular adsorption is found to be 0.05 eV and 
0.19 eV per water molecule more stable than partial or 
full dissociation, respectively. In view of the very contra- 
dictory results from previous DFT calculations in the lit- 
erature on the relative stability of molecular and dissocia- 
tive adsorption of water on TiO2(110), we have carefully 
re-evaluated the convergence of our computational setup. 
From the surveil of the literature it appears that the slab 
thickness and the degree of atomic relaxation of the slabs 
are the most crucial parameters. Therefore we have re- 
calculated the energy difference AE^d between molecular 
and dissociative adsorption of water as a function of slab 
thickness. In the first set of calculations we have consid- 
ered a full monolayer water coverage as represented by 
one water molecule in a (1x1) surface unit cell. As shown 
in Fig. [ini the molecular and dissociative adsorption en- 
ergy ^^'^ ^ad^' respectively, exhibit the well-known 
strong odd-even oscillations if the water molecules are 
adsorbed symmetrically on both sides of the slab and a 
full atomic relaxation is performed. These oscillations are 
even present in the energy difference AEe,d=E^^^—E^^^. 
Switching to our standard setup in which we adsorb wa- 
ter only on one side of the slab, fix the atoms in the two 
bottom layers at the bulk positions and saturate the bro- 
ken surface bonds at the bottom with pseudo-hydrogen 
atoms, the odd-even oscillations are almost completely 
removed (see Fig. [TU|) . For Ai?ad we find only a small 
increase from 0.19 eV for the four trilayer slab to 0.23 eV 
in the limit of infinite slab thickness. In the second set 
of calculations we repeated the calculations for a single 
water molecule in a (4x2) surface unit cell with our pas- 
sivated slabs for slab thicknesses up to seven trilayers. In 
the seven trilayer calculation the supercell contained 355 
atoms. Also in this case, as shown in Fig. [TOl there are 
almost no odd-even oscillations visible in the water ad- 
sorption energy i?ad and the energy difference Ai?ad be- 
tween the molecular and dissociative adsorption increases 
only slightly from — O.lleV for four trilayers to — 0.06eV 
in the six and seven trilayer calculation. Since in gen- 
eral such small energy differences depend sensitively on 
the quality of the pseudopotentials and the functional 
which are used in the computations, we only conclude 
that molecular and dissociative adsorption become ener- 
getically degenerate in the limit of single, isolated water 
molecules. 

When comparing our results on the adsorption of iso- 
lated water molecules with the literature surveil in Ta- 
ble |V] it has to be taken into account that in the previous 
calculations for 'single' water molecules smaller surface 
unit cells with effective water coverages of 1/4 monolayer 
were used, compared to 1/8 monolayer in our study. As 
can be seen from our results in Table IVIl at 1/4 mono- 
layer coverage there is already a significant interaction 
between the water molecules. For two water molecules 
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FIG. 10: Convergence of the molecular and dissociative water 
adsorption energy _Ead and the difference Ai5ad=i5S°' ^E^d^ 
with slab thickness at full monolayer coverage (left) and for 
single water molecules in a (4x2) surface unit cell (right). 
The solid black lines represent results from full relaxations 
with water adsorbed symmetrically on both sides of the slab. 
The calculations in which water was adsorbed only on one 
side of the slab with the bottom two trilayers fixed at the 
bulk positions and the broken surface bonds at the bottom of 
the slab saturated with pseudo-hydrogen atoms are shown by 
the red dashed lines. 



in the (4 x 2) surface unit cell the water binding energy is 
clearly reduced compared to a single molecule. The rea- 
son is, as we will point out below, the large contribution 
of the surface relaxations to the binding energy, which is 
now smaller since less degrees of freedom are associated 
with each molecule. In particular, for a (2x2) arrange- 
ment of water molecules, which was used in basically all 
previous calculations as setup for 'isolated' molecules, 
we find already a preference for molecular adsorption. 
Clearly, the 1/4 monolayer coverage does not represent 
isolated water molecules, and in principle, it would have 
to be checked if our 1/8 monolayer truly represents the 
limit of independent molecules. This implies that also in 
the embedded cluster approach large clusters are needed 
to take these contributions from the atomic relaxations 
appropriately into account. 

Our result that for isolated water molecules molecu- 
lar and dissociative adsorption is energetically degener- 
ate is not in agreement with the common view expressed 
in most experimental studies that water only dissociates 
at defect sites. Since our calculations are very well con- 
verged, this might point towards a deficiency of the pseu- 
dopotential approximation or, in general, PBE/DFT cal- 
culations. However, it also has to be taken into account 
that all experiments have shown a certain amount of dis- 
sociated molecules and it is very difficult to distinguish 
whether they only stem from dissociation at defects or if 
some initial dissociation at very low coverages has taken 
place. 

Comparing the water binding energies of two water 
molecules in a c(4x2) with the p(4xl) arrangement we 
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find that the interaction of water molecules between dif- 
ferent [001] rows is rather weak since the relative position 
of the water molecules within the [001] rows does not 
matter. This is very different from the situation when 
both water molecules are placed into the same [001] row 
(as in the p(2x2) configuration), where the changes in 
the surface relaxation lead to a strong repulsion between 
dissociated molecules. This repulsion has also been re- 
ported by Lindan et al.^ although with a slightly smaller 
magnitude. The repulsion can be compensated by ad- 
sorbing the two water molecules on neighboring Ti(5) 
sites, since now a hydrogen bond between the molecules 
is formed. We have considered all possible pairs of disso- 
ciated and undissociated water molecules with 'upward' 
and 'downward' orientation as illustrated in Fig. [9l As 
can be inferred from the binding energies in Table IVll all 
these configurations are almost equivalent, with a slight 
preference of pairs of dissociatcd/undissociated molecules 
in which the donating molecule of the hydrogen bond 
stays intact and the receiving molecule dissociates. Fur- 
thermore, we can see that the relative orientation of the 
molecules within the pair ('upward' or 'downward') does 
not matter, which will be also true at higher water cov- 
erages. 

We take now these water pairs as building blocks to 
form higher water coverages. First we put two water pairs 
in the same [001] row to built structures of alternating 
full and empty [001] lines. The molecular adsorption is 
now the most stable one. The half-dissociated structure 
of Lindan et al.'^^ is only slightly lower in energy. The 
fully dissociated structure, however, has become clearly 
unfavorable. Since the interaction of the water molecules 
between [001] rows is weak, as already pointed out above, 
this result remains unchanged when all [001] rows are 
filled with water molecules and the full monolayer cover- 
age is reached. 

For the length of the Ti(5)-0water bond we obtained 
values between 2.246 A and 2.306 A for isolated water 
molecules and at full monolayer coverage, respectively. In 
a recent quantitative analysis of photoelectron diffraction 
experiments^^ a value for this bond length of 2.21±0.02 A 
was determined. Taking into account that GGA typi- 
cally overestimates bond length by 1-2%, our result for 
isolated water molecules agrees quite well with the ex- 
perimental value, but overestimate it slightly in the limit 
of full monolayer coverage. 

In conclusion, we find that in the limit of very low cov- 
erages the molecular and dissociative adsorption of water 
is almost degenerate. For single, isolated water molecules 
we find a small preference for dissociation, however, al- 
ready at coverages when water pairs can form, mixed 
dissociated/molecular structures are more favorable than 
full dissociation. Finally, in the full monolayer coverage 
limit the molecular adsorption is the most stable config- 
uration. Based on these results one would expect that if 
water islands are formed upon adsorption, almost all wa- 
ter molecules stay intact and only those molecules at the 
boundary of the water patch which do not have a neigh- 



bor to which they can form a hydrogen bond may disso- 
ciate. It is interesting to see the difference in the role of 
donated and accepted hydrogen bonds for the stability of 
water molecules between the TiO2(110) and ZnO(lOlO) 
surface. On ZnO(lOiO), water molecules also show a ten- 
dency to dissociate But while on ZnO(lOlO) iso- 
lated water molecules stay intact and only start to disso- 
ciate when they receive an hydrogen bond, isolated water 
molecules on TiO2(110) are found to be close to disso- 
ciation and recombination is preferential when they can 
donate a hydrogen bond. 

Our results for the water binding energies are 
throughout slightly lower than the previous results ob- 
tained with the PW91 functional^iiiii^i^ but signif- 
icantly higher than what has been found in RPBE 
calculationali^i^iiiii^iSS (see also Tables |IV] and |V|. In 
order to compare our water binding energy of 0.82 eV 
at full monolayer coverage with experiment, we have fur- 
thermore calculated the corrections due to the zero-point 
vibration energy (ZPE). The vibration frequencies of the 
surface with and without the molecular water monolayer 
were obtained in harmonic approximation by a finite dif- 
ference scheme. Only the atoms in the adsorbate and 
the first surface layer were displaced by 0.01 A in the 
three cartesian directions. For the ZPE we obtain 0.12 eV 
so that we arrive at a ZPE-corrected water binding en- 
ergy at monolayer coverage of 0.70 eV, which is in excel- 
lent agreement with the results of the TDS experiments 
of 0.65 eV— and 0.66 eV--. The minor overbinding of 
the water molecules is typical for PBE calculations and 
has also been observed for the adsorption of water on 
ZnO(lOlO), where PBE calculations gave a water binding 
energy of 1.13 eV (without ZPE corrections) compared to 
the TDS/Redhead analysis of 1.03 eV^ 

Looking at the qualitative trend of the adsorption en- 
ergy with water coverage of the surface it is interesting 
to see in Table IVII that overall the binding energy de- 
creases with water saturation. This contradicts the re- 
sults of Ref. 46, whereas it is in full agreement with the 
RPBE calculations of Lindan et al^i^ and, in addition, 
to the experimental observation in the TDS measure- 
ments. This is a rather unusual behavior since one would 
expect an increase in the binding energy with water cov- 
erage due to the formation of hydrogen bonds. In order 
to gain deeper insights into the origin of this behavior 
we have decomposed the process of desorbing the water 
molecules into four steps. Starting from the fully relaxed 
adsorbate structure we detach in the first step the water 
layer from the surface without relaxing the atoms. Then 
we separate the water molecules in order to obtain the 
hydrogen bond strength. Finally, we allow the TiO2(110) 
surface and the water molecules to relax to obtain the re- 
laxation energies of the systems (see Ref.^^ for a similar 
analysis of the water adsorption energy on ZnO(lOlO) ). 
For structures with neighboring water molecules we find 
indeed that hydrogen bonds are formed which increases 
the water binding energy. This gain in adsorption en- 
ergy, however, is overcompensated at higher coverages 
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by a loss in relaxation energy of the substrate per adsor- 
bate molecule. Why are the surface relaxations so much 
more important for the TiO2(110) surface than, for ex- 
ample, ZnO(lOlO)? As we have seen in Sec. |TT1 the sur- 
face relaxation energy of the stoichiometric, ideal surface 
is quite large, 1.21 eV per surface unit cell compared to 
0.37 eV for ZnO(lOlO). Thus, re-relaxations induced by 
the restoration of the cation bulk coordination upon wa- 
ter adsorption will have a much stronger impact on the 
binding energy of the molecules than for ZnO(lOlO). 

Lindan et al^ have suggested that the dissociation of 
water at low coverages has not been observed experimen- 
tally because the dissociation of the molecules is hindered 
by a significant dissociation barrier which is difficult to 
overcome on the time scale of the experiments. Using 
NEB with a subsequent refinement of the transition state 
with the dimer method we have calculated the activation 
barrier for water dissociation for two cases: an isolated 
water molecule in a (4x2) surface cell and a single water 
molecule within a molecularly adsorbed water layer (also 
employing a (4x2) supercell). Our results of 0.16eV and 
0.14 eV, respectively, are significant lower than the ac- 
tivation barriers reported by Lindan et al^ of 0.45 eV 
(single water molecule in a (4x1) cell) and 0.26 eV (dis- 
sociation of every second water molecule in a full mono- 
layer) and show a much weaker dependence on the water 
coverage. Therefore, based on our results, we cannot 
confirm that the dissociation of water molecules at low 
coverages is prevented by a large dissociation barrier. 



E. Surface Phase Diagram 

Finally we combine the results of all previous sections 
to construct a two-dimensional phase diagram. We as- 
sume that the TiO2(110) surface is simultaneously in 
thermodynamic equilibrium with reservoirs with which it 
can exchange O and H atoms, for example, a surround- 
ing O2 and H2 gas phase. The surface free energy A7 of 
the different TiO2(110) surface structures now depends 
on both chemical potentials, A/iQ and A/ijj- By indi- 
cating the most stable structure and composition of the 
TiO2O(110) surface as a function of the two chemical po- 
tentials we obtain the phase diagram which is shown in 
Fig.ini 

The surface phase diagram of Fig. [11] summarizes in a 
condensed fashion the results of our study on the rela- 
tive stability of the different surface compositions which 
we have considered. The phase diagram is dominated by 
four different types of surface structures: the stoichiomet- 
ric, ideal surface, the water saturated, oxidized surface 
and the reduced surface structures by O depletion and H 
adsorption. For ambient conditions (A/^o = — l.leV) 
we find that the TiO2(110) surface is saturated with 
molecularly adsorbed water. The water can be gradu- 
ally removed by lowering the O and H chemical poten- 
tials (i.e., by heating or lowering the partial pressures). 
Passing through a transition region with low water cov- 
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FIG. 11: Phase diagram for the TiO2(110) surface in ther- 
modynamic equilibrium with H2 and O2 particle reservoirs 
controlling the chemical potentials A/ih and AfiQ. 



erages in which dissociated water molecules may appear, 
the stoichiometric, ideal surface structure is reached. 
It is the most stable surface structure at UHV condi- 
tions for a wide range of temperatures (A/iq = — 1.9eV, 
A/xh = — 1.5eV at room temperature). By heating the 
surface in UHV, the surface can be reduced by removing 
bridging oxygen atoms and creating structural O vacan- 
cies. At hydrogen rich and oxygen poor conditions the 
surface is reduced by hydrogen adsorption. However, the 
surface cannot be fully hydroxylated, but only a maxi- 
mum H coverage of about 60% can be reached. 

Overall the surface phase diagram shows a quite ex- 
pected behavior of the TiO2(110) surface and does not 
yield any new surprises. Under UHV conditions the sto- 
ichiometric and the O reduced surface structures are the 
most important ones. On the other hand, under wet 
conditions or at high hydrogen partial pressures, as it 
is typical in catalytic or photochemical applications, the 
surface is water covered or reduced via hydroxylation. 



IV. SUMMARY AND CONCLUSIONS 

Using density functional theory we have investigated in 
detail the formation of oxygen-related vacancies as well as 
the interaction of hydrogen and water with the TiO2(110) 
surface in terms of a comprehensive surface phase dia- 
gram. Careful convergence tests have shown that a com- 
putational setup of using slabs consisting of four trilay- 
ers in conjunction with saturating the surface bonds at 
the bottom with artificial hydrogen-like atoms with ef- 
fective nuclear charges of -1-4/3 and -1-2/3 and fixing the 
atoms in the bottom two trilayers at their bulk positions 
almost completely eliminates the otherwise pronounced 
odd-even oscillations of surface and adsorption energies 
as a function of slab thickness. Using this procedure, the 
relaxation energy of the stoichiometric surface as well as 
the H and H2 O adsorption energies at full monolayer cov- 
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erage are within 0.06 eV per surface unit cell, 0.02 eV per 
H atom and 0.04eV per H2O molecules, respectively, of 
the extrapolated results for infinitely thick slabs. Thus, 
very well converged results with slab thickness are ob- 
tained without relying on somewhat questionable extrap- 
olation schemes such as the "25% rule" 4- or averaging 
four and five trilayer results, '^^ as proposed previously in 
the literature. 

It is also demonstrated that the contribution of sur- 
face relaxations to both binding and formation energies of 
molecules and defects, respectively, are substantial. Both 
adsorbate- and defect-induced relaxations are not limited 
to one surface unit cell, but contributions to the relax- 
ation energy from more neighboring surface cells have to 
be taken into account. This leads to substrate- mediated 
interactions between defects and adsorbate molecules 
which can extend over several lattice constants. Con- 
sequently, large surface unit cells have to be used when 
properties of isolated defects or adsorbed molecules are 
studied. The same arguments apply for embedded clus- 
ter studies where clusters have to be large enough to be 
able to take the relaxations of the relevant environment 
appropriately into account. 

In general, GGA/DFT calculations of reduced Ti02 
structures are hampered by the underestimation of the 
Ti02 band gap of about 1.0 eV. For isolated O vacan- 
cies and pairs of OH groups on the TiO2(110) surface it 
has been shown that in GGA/DFT calculations the de- 
fect level lies at the bottom of the conduction band^i^ 
instead of 0.7-0.9 eV below. Naively one would there- 
fore expect that GGA/DFT calculations should underes- 
timate the energy of reduced Ti02 structures and thus 
overestimate the reducibility of Ti02 since defect states 
will tend to be too delocalized. In contrast to this ex- 
pectation, a survey of the literature data showed instead 
that GGA/DFT overestimates the O vacancy formation 
energy. Furthermore we found in the present study that 
the dissociation energy of water at O vacancies is slightly 
too low in the PBE/DFT calculations and also the pre- 
dicted maximum coverage of the TiO2(110) surface with 
hydrogen atoms is found to be slightly smaller than ob- 
served in experiment. This implies that in all three cases 
the GGA/DFT calculation have underestimated the re- 
ducibility of the TiO2(110) surface. However, the pre- 
sumable error in the total energies is small: about 0.3 eV 
for the O defect formation energy and a tenth of an eV 
in the other two cases. Based on this analysis we con- 
clude that the energetics is quite reliably described by 
GGA/DFT despite the problem of such calculations to 
obtain accurately the underlying electronic structure of 
the corresponding reduced Ti02 surfaces. It might be ar- 
gued that the energetics of reduced Ti02 is well described 
because although the defect level of O vacancies and OH 
pairs is at the bottom of the conduction band instead 
of 0.7-0.9 eV below, its energy difference with respect to 
the top of the valence band is almost right. This might 
be correct for single, isolated O vacancies and OH pairs. 
However, in the case of an increasing hydroxylation of the 



surface when the defect states start to interact and to hy- 
bridize with the conduction band the whole defect band 
is now more easily accessible due to the underestimation 
of the band gap, but still the overall energetics is quite 
reliable with a tendency to underestimate the reducibility 
of Ti02. 

After having clarified these crucial convergence and ac- 
curacy issues, we have compared the relative stability 
of surfaces with oxygen vacancies or adsorbed hydrogen 
atoms and water molecules by combining our DFT calcu- 
lations with a thermodynamic formalism. Lowest-energy 
structures were determined by assuming that the sur- 
faces are in thermodynamic equilibrium with oxygen and 
hydrogen reservoirs. The exchange of O and H atoms 
with the reservoirs is described by introducing appropri- 
ate chemical potentials which can be related to exper- 
imental temperature and pressure conditions. This al- 
lows us to extend the zero-temperature and zero-pressure 
DFT calculations to more realistic environmental condi- 
tions such as in UHV surface science experiments or in 
catalytic and photochemical applications. 

The central result of this approach is the phase di- 
agram depicted in Fig. 1111 For ambient conditions we 
find that the non-reductive adsorption of water prevails. 
A full monolayer of molecular adsorbed water is pre- 
dicted as the most stable adsorbate structure in agree- 
ment with experimental observations. In the low cover- 
age limit molecular and dissociative adsorption becomes 
energetically degenerate. In addition, according to our 
calculations the dissociation barrier for molecularly ad- 
sorbed water molecules is not large enough to prevent 
their dissociation even at low temperatures. 

At strongly reducing, oxygen-poor conditions struc- 
tural O vacancies are formed on the TiO2(110) surface in 
thermodynamic equilibrium. However, for a wide range 
of temperature and pressure conditions it is found that 
removing full Ti02 units is thermodynamically more fa- 
vorable instead. In the case of reduction of the surface via 
hydrogen adsorption we observed that the surface cannot 
be fully hydroxylated in thermodynamic equilibrium. A 
maximum coverage of about 60% is predicted from our 
thermodynamical analysis. Furthermore, our calculated 
energy barriers for hydrogen migration indicate than hy- 
drogen atoms rather migrate into the bulk than desorb- 
ing from the surface. Both observations are supported 
by recent experiments in which a maximum saturation 
of the surface with hydrogen of 70% is reported, whereas 
no desorption of H2 could be observed in TDSi^ 
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